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Abstract. The aim of this article is to promote the use of probabilistic methods in the 
study of problems in mathematical general relativity. Two new and simple singularity 
theorems, whose features are different from the classical singularity theorems, are proved 
using probabilistic methods. Under some energy conditions, and without any causal or 
initial/boundary assumption, simple conditions on the energy flow imply probabilistic in- 
completeness. 



The introduction of the concepts of singularities and boundary of a spacetime both have 
their source in the observation that Einstein's geometrical picture of a spacetime does not 
prevent the existence of some undesirable features like the existence of regions where some ge- 
ometrical scalar explodes along some path, or the existence of physically relevant incomplete 
paths (geodesies, paths with bounded acceleration, etc). Following the pioneering works of 
Penrose [T| and Hawking [2], [3], [I], most singularity theorems state that a spacetime has 
an incomplete causal geodesic provided some energy, causality and boundary or initial con- 
ditions hold. We consider in this article a purely geometrical random dynamics producing 
random timelike paths, to be considered as randomly perturbated geodesies. Technically 
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speaking, we are going to construct some non-trivial probabilities on the (separable metric) 
space of inextendible timelike paths parametrized by their proper time, with finite or infinite 
lifetime. We prove in section [3] that some energy conditions and conditions on the energy 
flow are sufficient to ensure the existence of a set of inextendible incomplete timelike paths of 
positive probability. This provides an unusual type of conclusion under a non-common set of 
conditions where no causality assumption is required, nor any boundary or initial condition. 

Cartan's moving frame language provides a very intuitive way of describing C 2 trajectories 
by transporting parallelly an initial frame along the path and describing the variations of the 
velocity in that moving frame. The datum of these vector space-valued variations suffices 
to reconstruct the path. Taking them as the primary object and choosing them randomly 
produces random paths (forgetting all possible technical problems). Our construction of 
random timelike paths relies on a variant of this procedure and associates to any starting 
point (mo, mo) i n the future-oriented unit bundle T^M a timelike path in M under the form 
of a continuous path (m s ,m s ) in T X M subject to the condition -^m s = m s . So our random 
dynamics is a dynamics in T^M, in the same way as the (timelike) geodesic motion is better 
viewed as a flow in the unit bundle. The (starting point dependent) distribution of these 
random timelike paths will provide the above mentionned set of probabilities. 

These random dynamics are called "diffusions" in the sequel. The reader should have in 
mind that these diffusionss are not physical diffusions in the sense that they are not models 
for diffusion phenomena in a relativistic medium like a gaz. We refer for this kind of questions 
to the works by Debbasch and his co-workers, [5], [6], or Dunkel and Hanggi, [7], [8], their 
review [9], and the references cited therein. Rather, the diffusions considered below are what 
probabilists call a diffusion: intuitively, a random process X = (X 1 , . . . ,X n ) such that 

E[Xi +e -Xi] =b\X s )e + o{e) ) 

E[(Xl +e - Xl)(Xl +e - X*)] = a*(X,)e + o(e). 

for i,j = l..n,s ^ and e > 0. Here, T s = c(X r ; r ^ s) and b l and a* J are measurable 
functions; see chapter 5 in |10| . 

Section [2] gives a brief overview of the diffusions considered here, with the help of which 
we derive two simple probabilistic singularity theorems in section |3l theorems [8] and [T2j As 
mentionned above, the random dynamics considered in this work are random perturbations 
of the geodesic flow in T^M. The main issue addressed in this work is the following question: 
Can (timelike) geodesic completeness be destroyed by a random perturbation of its dynamics? 
Theorems [8] and [12] both provide conditions under which the answer to this question is 
positive. The core of their proofs is a stochastic analogue of the following kind of trivial 
observation about the unperturbed flow. If there exists two functions / ^ h and an initial 
condition <po G T^M such that / is bounded below by e cs and h bounded above by e c s along 
the geodesic started from </>o, then this geodesic cannot be future complete if d < c. The 
stochastic counterpart of this observation is subtler though. 

As this work is the first ever written on the subject, we have chosen to present only some 
basic aspects of the situation. We hope this will convince the reader of the possible intest 
of using probabilistic methods in problems about spacetime geometry. 
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2. Relativistic diffusions 

One owes to Dudley [TT] and Schay |12| . |13| . the merit to have asked first if there is a 
natural way of defining random timelike paths in Minkowski spacetime, without introducing 
any additional structure than the metric. The models of physical diffusions considered e.g. by 
Debbasch et al. j5], Angst and Franchi |14] , or Dunkel and Hanggi |7], require an additional 
vector field to be defined. Dudley provided a complete answer to the above problem but 
could not pursue further his investigations to the general relativistic framework for lack of 
technical tools available at that time. The subject remained untouched for nearly forty years 
before Franchi and Le Jan generalized Dudley's motion to a general Lorentzian manifold in 
|15j . We give a brief overview of these dynamics in section 12711 and [2721 below as they will be 
our main tool to probe the existence of spacetime probabilistic singularities in section [3j 

2.1. Relativistic diffusion in Minkowski spacetime. The question asked by Dudley 
and Schay is the following: Is there a natural way of defining Markovian random time- 
like paths, independently of any reference frame, and without using any other object than 
Minkowski metric? To describe their answer, denote by g Minkowski metric and write 
M= {me E 1 ' 3 ; m° > 0, g(m,m) = l} for the upper half-sphere. Although g has signature 
(— , +, +, +), its restriction to any tangent space of the spacelike hypersurface EI is definite- 
positive, so HI inherits from the ambient space a Riemannian structure which actually turns 
it into a model of the 3-dimensional hyperbolic space. Brownian motion on H is defined as 
the unique continuous Markov process {m s } s ^o with generator half of the Laplacian A of H. 
(We have indexed Brownian motion by E + ; it can indeed be proved that it has an infinite 
lifetime; see e.g. |16j.) As any timelike path {m s } s ^o indexed by its proper time is differen- 
tiable at almost-all times, and determined by its derivative m s , since m s = itiq + J^m u du, 
a random timelike path is determined by an H-valued random process {m s } s ^Q. So, if one 
wants to talk of a random timelike paths-valued Markov process, one needs to record the 
position and the velocity of the path in the state space and work in E 1 - 3 x H. Dudley showed 
in [11] that there exists essentially a unique E 1 ' 3 x H-valued continuous Markov process 
{(m s ,m s )} s ^ such that 

• one has m s = m + J Q S m u du, for all s ^ 0, 

• its law is invariant by the action of the affine isometries of E 1 ' 3 . 

It corresponds to a Brownian velocity process {m s } s ^>o with generator ^- A, where a is a 
positive constant. In short, there is (essentially) a unique way of constructing a random 
timelike path, by imposing to the velocity to undergo Brownian oscillations in EI. 

It will clarify the construction in the general framework of a Lorentzian manifold to give 
a slightly different picture of Dudley's process. Consider the random motion in spacetime 
of an infinitesimal rigid object, represented by a path in E 1 ' 3 x 500(1, 3), where we write 
6*00(1, 3) for the identity's component in 5*0(1, 3). Denoting by (m, e) = (m, (eo, e,\, e2, 63)) 
a generic element of E 1 ' 3 x 50o(l,3), the map 7Ti((m,e)) = (m, eo) is a projection from 
E 1 ' 3 x 5O (l, 3) onto E 1 ' 3 x EI. Denote by {e , e±, e 2 , e 3 } the canonical basis of E 1,3 . For each 
j G {1, 2, 3}, the Lie element Ej = eo ® e* + 6j <S> Cq e so(l, 3) generates a hyperbolic rotation 
in the 2-plane spanned by e and ej. Define on 50o(l, 3) some left invariant vector fields In- 
setting 

V j (e) = eE j , j €{1,2,3}. 
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Let w be a 3-dimensional Brownian motion and use the notation od for Stratonovich differ- 
ential. Consult section 5, chap. V, of the book |10J . or chapter V of |17j . for the necessity 
of using Stratonovich formalism when dealing with stochastic differential equations on man- 
ifolds, and for a comparison with Ito's differential. By construction, the solution (e s ) s ^o of 
the Stratonovich stochastic differential equation on 5*00(1,3) 

3 

(2.1) ode s = ^V j {e s )odw j s 

i=i 

projects down by 7Ti into a Brownian motion on H. (Consult the book [18] of Elworthy, or 
|17] , if you do not feel comfortable with this fact.) So Dudley's dynamics is described by 
equation (12. ip and 



(2.2) dm s = e (s) ds. 

Remark. It is interesting to notice that Dudley's diffusion is obtained heuristically as the 
large scale limit of geodesies in the discrete random approximations of a spacetime introduced 
by Sorkin in his causal set theory; see [19], [20] and |21| . 



2.2. Basic relativistic diffusion on a general Lorentzian manifold. Let (M, g) be a 

Lorentzian manifold, oriented and time-oriented. (These assumptions are harmless as they 
hold on a finite covering of M.) Denote by X^M the unit future-oriented bundle over M 
and by OM the component of the orthonormal frame bundle made up of pairs (m, e), with 
m G M and e = (eo, ei, e2, 63) a direct orthonormal basis of T m M with eo future-oriented. 
We shall denote by <fi = (m, m) a generic point of T X M and by $ = (m, e) a generic element 
of OM. Write m : OM T X M and ir : OM M for the canonical projections. 

As can be expected from the Minkowskian picture, the basic relativistic diffusion on a 
Lorentzian manifold M is actually a diffusion in T X M. Roughly speaking, one can con- 
struct this process by rolling without splitting the trajectories of the relativistic diffusion 
in Minkowski spacetime on M. A more formal way of proceeding is to introduce an OM- 
valued S'0(3)-invariant diffusion process whose projection in T X M is consequently a diffusion 
process on its own. 

The choice of OM as a framework is motivated by the fact that it bears more structure 
than T^M and is the natural framework where to use Cartan's ideas on moving frames. 
The action of S*Oo(l,3) on each fiber of 7r induces the canonical vertical vector fields; 
denote by Vj the vector field associated with the Lie element Ej. Denote by (-ffi)j=o..3 the 
canonical horizontal vector fields on OM associated with Levi-Civita connection. The basic 
relativistic diffusion is defined in a dynamical way as the unique solution of the following 
Stratonovich stochastic differential equation on OM 

3 

(2.3) orf$ s = H ($ s ) ds + a Vj ($ fl ) odw 3 ., 

where w is a 3-dimensional Brownian motion and a a positive constant. It describes a random 
perturbation of the geodesic flow whose intuitive meaning is the following. To get & s +ds out 
of $ s = (m s ,e s ), transport first e s parallelly along the geodesic starting from m s in the 
direction eo(s), during an amount of time ds; you get an orthonormal frame f s of T ms+ds M. 
Making then, in each spacelike 2-plane of f s , independent hyperbolic rotations of angle a 
centered normal random variable with variance a 2 ds, you get G s +ds- This dynamics is the 
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straightforward generalization of Dudley's dynamics (12. ip and (12. 2p and was first considered 
in [15]. The following statement gives a different view on this dynamics. 

Lemma 1. Let 7 : [0,T] — > M be a C 2 timelike path parametrized by its proper time, and 
T e OM such that ni(T ) = (7(0), 7(0)) G T X M. Then there exists a unique C 2 path 
{^s) 0<s<T i n OM, and some unique C 1 real-valued controls h l ,h 2 ,h 3 defined on [0, T], such 
that i = T , 7Ti(* s ) = (7(5), 7(s)) and 

3 

So the relativistic diffusion dynamics is obtained by replacing the deterministic controls of 
a typical C 2 timelike path by Brownian controls. 

Equation (12. 3p has a unique strong solution defined up to its explosion time (. Write OM 
as an increasing union of relatively compact open sets O n and define the stopping times 
T n = inf{s ^ ; <f> s ^ O n }. The explosion time ( is by definition the increasing limit of the 
T n , and does not depend on the arbitrary choice of sets O n . It is not difficult to see on this 
equation that it defines an S'Oo(3)-invariant diffusion which, as a consequence, projects on 
T 1 M into a diffusion process. Consult |15| . theorem 1, or [22], theorem 4, for the details. 
The diffusion (^s) 0<s< ^ i n OM has generator 

where we consider vector fields as first order differential operators. In so far as we are 
primarily interested in the T 1 M-valued process (as it provides us directly with random 
timelike paths), while we shall mainly be working with the above OM- valued process it is 
important to notice that 

PROPOSITION 2 ([23], prop. 1). The OM-valued diffusion ($s) 0<s<? and its m projection 
in T : M have the same lifetime. 

We shall thus freely work in the sequel with the OM-valued diffusion {& s )o^s<(- 

Remarks, (i) You might ask why we called basic relativistic diffusion the solution of 
equation (12. 3p . and not simply relativistic diffusion. This is due to the fact that, contrary 
to what happens in Minkowski spacetime, where there is only one way of constructing nice 
random timelike paths, there are many ways of doing it on a general Lorentzian manifold. 
Think for example of a diffusivity a in (12. 3p depending on the location of the particle (it 
may be the scalar curvature of the manifold at that point for instance). Consult [22J and 
|24j for more material on this subject, seen from a mathematical point of view, and [25] for 
a physical point of view on related matters. Let us repeat here that this diffusion is not to 
be thought of as a mathematical model for a physical diffusion phenomenon but rather as a 
mathematical object useful for studying some features of the spacetime geometry. 

(ii) Basic relativistic diffusions have only been studied explicitl in a few cases: in Minkowski 
spacetime [2S], Robert son- Walker spacetimes [27J, Schwarzschild [T5] and Godel [2B] space- 
times. The stochastic completeness question is easily delt in each case. While it is trivial to 
find a geodesically incomplete spacetime which is stochastically complete (remove a point 
from Minkowski spacetime - see paragraph 13.11 below), the possibility to have a (timelike) 
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geodesically complete and stochastically incomplete spacetime has noot been established so 
far; theorems [8] and [12] address that issue. 

(iii) As mentionned in the introduction, one can consider equation ( 12. 3ft as a dynamical 
way of constructing a probability measure on the separable metric space of inextendible time- 
like paths parametrized by their proper time (with finite or infinite lifetime). Given a starting 
point $ £ OM, the distribution P$ of ($ s )o^s<c i s such a probability. Each such probabil- 
ity is highly non-trivial, as shown by the following qualitative statement. Let 7 : [0, 7} — > M 
be a C 2 timelike path parametrized by its proper time, and (ei(0), e2(0), 63(0)) be such that 
(7(0),ei(0),e 2 (0),e 3 (0)) is a frame of T 7(0) M. Transporting parallelly (ei(0), e 2 (0), e 3 (0)) 

along 7 provides a frame (7(5), ei(s), e 2 (s), e 3 (s)) of T 7 ( S )M. The map 

3 

F : (s,x) E [0,7] xl 3 4 exp^^Ve^s)) 

i=i 

is a well-defined diffeomorphism from [0,7] x U onto its image V, for some small enough 
open ball U of M 3 with center 0; V is a tube around 7. 

PROPOSITION 3. Let $ = (m, e) e OM be such that 7(0) e I + (m). 7hen the IF '^-probability 
that the basic relativistic diffusion hits -F({0} x UJ and exits the tube V in F({T} x U\ is 
positive. 

An example: basic relativistic diffusion in Schwarzschild spacetime. Let denote by 
(M, g) Kruskal-Szekeres extension of Schwarzschlid spacetime. Franchi and Le Jan proved 
in their seminal work [15J that, for any starting point, the relativistic diffusion hits the 
boundary of the black hole with a positive probability strictly smaller than 1. After that 
time, as any other timelike path, the trajectory of the diffusion hits the singularity before 
units of proper time have ellapsed, where R is the radius of the black hole. So ( is finite 
with positive probability, whatever the starting point of the diffusion. 

3. Probabilistic incompleteness 

3.1. Probabilistic incompleteness and geometry. A probabilistic incompleteness the- 
orem is a statement of the form P$(C < 00) > 0, for some initial starting point $ G OM of 
the diffusion process. It should be clear from the description of Dudley's diffusion given in 
section 12.11 that probabilistic incompleteness and geodesic incompleteness are two different 
notions in general. Indeed, removing a point from Minkowski spacetime produces a geodesi- 
cally incomplete spacetime while Dudley's diffusion is almost-surely defined for all proper 
times. This is easily seen as follows. Choose a frame and denote by t s the time component 
of the M 1,3 -part of the diffusion. As -£-t s ^ 1, one can change parameter and use the time 
of the frame rather than the proper time to parametrize the trajectory. This amounts to 
consider the diffusion as an (JBL 3 x H) -valued diffusion. It is easily seen to be hypoelliptic, 
using Hormander's criterion, so that its distribution at any given time has a density. If the 
point removed belongs to the slice {t = 7}, the diffusion has a null probability of hitting it. 
It is also true that any 2-dimensional submanifold of OM is polar for the basic relativistic 
diffusions. We describe in theorems [8] and [12] below geodesically complete spaces which are 
probabilistically incomplete. 

In a Riemannian setting, it is well-understood that any statement about the Laplacian, 
its eigenvalues/eigenfunctions, heat kernel, zeta function, etc. provides a non-elementary 
information about the geometry of the manifold. In so far as the generator Q of the basic 
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relativistic diffusion is constructed out of canonical geometric vector fields, any statement 
about Q provides a non-elementary information about the geometry of OM, or T X M. 

The following proposition links probabilistic incompleteness and geometry; its proof is 
similar to the proof of the corresponding statement for the Laplacian, and can be found e.g. 
in theorem 6.2. of [29]. 

PROPOSITION 4. The following statements are equivalent. 

(1) There exists a point $ G OM such that P$(£ < oo) > 0. 

(2) Let A > 0. There exists a non-null bounded function f : OM — > R such that 



(3) Let T > 0. There exists a non-null bounded solution to the Cauchy problem 
dth = Qh, on [0,T] x OM, with initial condition 0. 

This equivalence between the probabilistic incompleteness problem and the above two prob- 
lems on linear partial differential equations brings a different point of view on the probabilis- 
tic explosion problem, and a huge tool kit to investigate it. There is no similar correspondence 
for the classical geodesic incompleteness problem. The following proposition clarifies what 
happens in case of explosion. 

PROPOSITION 5. Suppose there exists a $ G OM such that P$(£ < oo) > O.Then explosion 
occurs in an arbitrarily small time with positive probability: P$(C < e) > 0, for all e > 0. 

PROOF - Set f(^f,s) = P^(C < s), for \1> G OM and s > 0; due to Hormander's theorem, 
this function depends smoothly on s) G OM x . Write $ = (m, e) and set 
/+($) = 7T- 1 (/ + (m)). Given e > 0, suppose P^(C < e) = for all $ G /+($). Using 
the strong Markov property first, then the Markov property inductively and the fact 
that C for G /+($), we would have P$(C > ke) = 1, for all k ^ 1, 

contradicting P$(C < oo) > 0. So there exists an element \l/o £ -^ + (^ ) ) such that 
/(^Oi e ) > ; by continuity, this remains true in an open neighbourhood of \l/o- As, 
by proposition [31 one can reach this open set in an arbitrarily small time from $ with 
positive probability, the result follows from the strong Markov property. > 

Probabilistic incompleteness results have two other noticeable distinct features when com- 
pared with the classical deterministic singularity theorems. As a function of $, the quantity 
P$(C < oo) is known to be (/-harmonic (Gf = 0) and satisfies the strong maximum principle. 
So if P$(C < oo) = 1 at some point $ = (m, e) G OM, we actually have P$'(C < oo) = 1 for 
any $' = (m', e') with m! in the chronological future of m. 

Given a spacelike hypersurface V of M denote by ay the volume measure on V inherited 
from the ambient geometry and set OV = = (m, e) G OM ; m G V}. Write Vol©v 
for the natural volume measure VoL m (c?e) ® ay {dm) on OV, where VoL m (<ie) is the Haar 
measure on the fiber SOq(1, 3) above m. Given a point $o = ( m (b e o) £ OM, there exists a 
relatively compact neighbourhood U of mo in M such that U is included in a larger globally 
hyperbolic neighbourhood of mo and d(lA H J + (mo)) is the union of a lightlike hypersurface 
and a smooth spacelike hypersurface V. Denote by H the random hitting time of OV by 
($ s )o<; s <(; it is P$ -almost-surely finite. The hypoellipticity of the generator Q of the diffusion 
garantees that the distribution of $r has a smooth positive density D(ty) with respect to 
Volqv. Consult proposition 5 of [22] for a proof. We can thus write 



(0-A)/ = O. 
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So if < P$(C < oo), there is a subset £ of OV of positive VoLov- meas ure such that 
P*(C < oo) > for all \P £ £. Not only do we have a set of positive probability of 
incomplete inextendible timelike paths with a common starting point, but this conclusion 
also holds for a non-trivial uncountable collection of starting points. 

3.2. Completeness conditions. Despite all the works done, it remains unclear what pre- 
cise features of a spacetime forbid the existence of honest inextendible incomplete timelike 
paths. The example found by Geroch (30) of a timelike and lightlike geodesically complete 
spacetime having an inextendible incomplete path with bounded acceleration gives an idea 
of the subtleties involved in that matter (see also [31]). 

As far as our method is concerned, the recent work [23] delimitates its domain of ap- 
plication by determining some general situations in which the paths of the diffusion have 
almost-surely an infinite lifetime. We recall here two such conditions and refer the reader to 
that work for more and subtler material. 

• Let M = I x S be a globally hyperbolic spacetimes whose metric tensor is of the form 

9m(q, q) = -a 2 m \q°\ 2 + h m (q s , q s ), q £ T m M, 

where q° is the image of q by the differential of the first projection / x S — > I and q the 
image of q by the differential of the second projection / x S — > I. Write m — (t,x) £ / x S. 
The function a is a positive C 1 function and h m a positive-definite scalar product on T X S, 
depending in a C 1 way on m. This class of spacetimes contains as elements all Robertson- 
Walker spacetimes - in particular de Sitter and Einstein-de Sitter spacetimes - and the 
universal covering of the anti-de Sitter spacetime. 

PROPOSITION 6 ([23]). The relativistic diffusion does not explode ifVa is everywhere non- 
spacelike, future- oriented. 

• It has been proved in [23] that if the spacetime is b-complete in the sense of Schmidt (a 
strong requirement; see e.g. [32] and [33]) then the relativistic diffusion does not explode. 

3.3. A first probabilistic singularity theorem, a) What can we expect? Pursuing 
the seminal works of Penrose and Hawking, the proofs of most singularity theorems follow 
one of the following two lines of reasonning (consult Senovilla's review [34] for a thorough 
and critical review on singularity theorems, or Wald's book |35|). 

(1) Using causality conditions, one constructs a useful maximal geodesic. At the same 
time, the energy and initial/boundary conditions induce conjugate or focal points 
along any geodesic, provided time is allowed to run long enough; as this cannot 
happen along a maximal geodesic, its time parameter has to be bounded. 

(2) Supposing the spacetime geodesically complete, one constructs a compact proper 
achronal boundary whose existence prevents the existence of an open Cauchy hyper- 
surface. 

Each of these schemes uses in a crucial way the rich structure of the geodesic flow. In so 
far as the random flow generated by equation (12. 3 P is not easy to analyse, it is not obvious 
to see which geometric features of a spacetime can lead to exploding solutions of equation 

(Q. 

The following very general and rough explosion result will be sufficient for our purpose in 
the next section; it is due to Khasminski. A generalisation of this result is proved in lemma 
[TT1 Recall the equivalence between points (1) and (2) in proposition [H 
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LEMMA 7. Let f : OM — > M. be a C 2 function bounded above by a positive constant and $ a 
point o/OM such that /($o) > 0- U there exists a positive constant C such that Qf ^ Cf, 
then P$ (( < oo) > 0. 

b) A probabilistic incompleteness theorem. Let R(-,-)- be Riemann curvature tensor 
and denote by RlC the Ricci curvature tensor and by RlC its restriction to T X M : Ric(0) : = 
RlC m (rh, rh), for = (m, m) E T X M. We identify Ric to a function on OM setting Ric($) = 
RlC (jti ($)) . Write R for the scalar curvature and T = RlC — ^ R g for the energy-momentum 
tensor. 

THEOREM 8. Let (M, g) be a Lorentzian manifold. Suppose 

(1) V (m, m) G T^M, T m (m,m)^0, 

(2) t/ie function RlC is non-identically constant, positive at some point $o, bounded 
above, 

(3) i/iere exists a constant c6l sitc/i i/iai ifo RlC ^ cRlC. 

Write c = C — 2 a 2 /or some positive constants C and a. Then the basic relativistic diffusion 
with diffusivity a 2 explodes with ¥^ -positive probability. 

The proof makes a crucial use of the following simple lemma. 

Lemma 9. We have Q RlC = H RlC + 2 a 2 RlC + 2 a 2 T. 

Remarks, (i) It is worth noting that theorem [8] holds regardless of any causality assump- 
tion. Hypotheses (1) and (2) are pointwise energy conditions and condition (3) a dynamic 
condition on the energy flow. This set of conditions is of a very different nature from the 
usual conditions of the classical singularity theorems. 

(ii) One cannot impose to RlC to be bounded, as Dajczer and Nomizu have proved in [36] 
that the spacetime is Einstein under this condition, so RlC is constant. The space S 2 x R, 
where S 2 is the 2-dimensional Lorentz manifold of constant sectional curvature 1, satisfies 
condition (2) and is not an Einstein manifold. 

(iii) The boundedness condition (2) seems too demanding for theorem [8] to be of any 
physical interest. Our aim here is more to promote a method by a simple example rather 
than by a technical work which would improve these conditions. We shall nonetheless see in 
section 13.41 how to weaken this assumption. The main probabilistic ingredient in the proof 
of theorem |S] is the basic explosion criterion stated in lemma [71 This is a very general and 
rough tool which does not take into account any peculiar feature of our problem. 

(iv) Let A be any subset of M. The past domain of dependence of A is the set D~(A) of 
points m of M from which any future directed timelike path starting from m eventually hits 
A. The future domain of dependence D + (A) is defined similarly using past directed timelike 
paths. The domain of dependence of A is D(A) = D~(A) UD + (A); it is globally hyperbolic. 
Note that if A is a relatively compact spacelike hypersurface, and m a point of D~(A), then 
there exists a constant T(m), depending on m, such that all timelike path starting from m 
hits A before proper time T(m). 

Given A C M, it is well-known that if the dominant energy condition holds on its domain 
of dependence D(A) and the energy momentum tensor vanishes on A then it vanishes on the 
whole of D(A), |37| . Lemma |9] provides for free a similar result under different hypotheses. 
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PROPOSITION 10. Let A be a relatively compact spacelike hypersurface of a spacetime (M,g). 
Suppose the following conditions hold on D~(A): 

• strong and weak energy conditions, 

• there exists a constant c G M such that H$ RlC ^ cRlC. 
If RlC = on A then RlC = on D~{A). 

Note that the strong and weak energy conditions together do not imply the dominant 
energy condition. 

PROOF - Combining the hypotheses of the proposition and lemma [9] we get the inequality 

G RlC ^ (a 2 + c) RlC ^ 0, 

for a big enough a. The function RlC is thus a non-negative ^-sub-harmonic function on 
OM. Given $ = (m, e) with m G D~{A), the stopping time H = inf{s ^ ; 7r($ s ) G A} 
is almost-surely bounded by a constant depending on $ , an d 

E$ [Ric($ H )] = Ric($) +EJ [ (0Ric) ($ P ) dr\ ^ Ric($) ^ 0, 

Uo J 

by optional stopping. As Ric($//) = this proves that Ric($) = and implies the 
result as $ is any point of tt- x (D~ (A)). > 

A similar result holds for D + (A) if HqKiC ^ cRlC for some constant cel. 
3.4. A second probabilistic singularity theorem. We keep the notations of section l3~3l 

a) Another explosion criterion. We are going to use in this section a refined version of 
lemma [7] for which no boundedness hypothesis is needed. 

LEMMA 11. Suppose there exists two non-null, non-negative smooth functions on OM, with 
f $C h, and two constants ^ d < c such that 

Qf ^ c f an d Gh ^ c' h. 

Let $o £ OWl be such that /($o) > 0. Then the basic relativistic diffusion started from $o 
explodes with positive probability. 

PROOF - Suppose C, is P$ -almost-surely infinite, so we can apply Ito's formula and write 
for any s > 

/($,) = f(%) + Mi + / gf($ u )du, 

Jo 

where is a P(j, () -martingale with respect to the filtration generated by the diffusion 
process. Taking expectation and applying Fubini, we thus have 

Eo> [/($,)] = /($ ) + f [Gf ($„)] ^ /($ ) + c f E $0 [/($„)] du, 

Jo Jo 

so Gronwall's lemma gives 

E* [/($.)] ^/We" 

Similarly, we have E$ [/i($ s )] ^ /i($o)e c s . We get a contradiction noting that, since 
/ ^ h, we should have 

/($ )e cs ^ E* [/($,)] < E* [/i($ a )] < M$o)e c ' s 
for all times, which cannot happen. > 
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We are going to apply this explosion criterion to the function / = RlC and a function h 
of the form RlC + U. We construct this function U in the next paragraph. 

b) Green function on EL Let denote by A the Laplacian on H and by G : EE x EI — > R+ the 
Green function of \ A. Due to the highly homogenous character of H, the quantity G(x, y) is 
actually a function of the hyperbolic distance from x to y. Let x re f G EI be any reference point 
and let p denote the hyperbolic distance function to x re f. A continuous function / : M — > R 
is said to be exp-bounded if sup x e~ ap ^ < oo for every a > 0. This definition does 
not depend on the choice of reference point x re f. It is well-known (consult for instance the 
classic book |38| of A. Friedman) that if / is exp-bounded the equation 

-Aw = -/ 
2 J 

has a unique solution which is null at infinity; it is given by the formula 

u(x) = J G(x,y)f(y) dy, 

where we write dy for the Riemann volume form of EL 

Let us come back on OM and identify, for each m G M, the future unit tangent bundle 
T^M to H by arbitrarily identifying an element of O m M to the canonical basis of R 1 ' 3 . The 

exp-boundedness character of the function RlC does not depend on this arbitrary choice; 
suppose it is exp-bounded for all m G M and set, for $ = (m, (e , • • • , e 3 )) 

(3.1) U(Q) = 2[ G(e ,y)RiC m (y,y)dy. 

Jt^m 

As G(eo, y) depends only on the hyperbolic distance from eo to y the function U is well- 
defined independently of our arbitrary identifications. It solves the equation 

i 3 

(3.2) -J2V 3 2 U = -2RlC, 

3=1 

3 

since the operator | Vf on induces on each fiber T^M the operator |A. 

3=1 

c) A second probabilistic singularity theorem. The following singularity theorem is 
similar in nature to theorem [HI and essentially states that a spacetime has a probabilistic 
singularity if some static and dynamical energy conditions hold. No causality assumption is 
needed. Set h = RlC + U, and recall we write R for the scalar curvature. 

THEOREM 12. Let (M.,g) be a Lorentzian manifold satisfying the following conditions. 

(1') Static energy conditions. RlC is non-negative and non-identically null, and R ^ 0. 
(2') Regularity condition. The function RlC is exp-bounded in eachT^M. and there exists 

some constants < a < 1, 0^c'<c and < a < \ such that 

V Z y 2 a 7 

1 — a 



Ric < U. 



a 



(3') Dynamical energy conditions, (i) H RlC ^ (c — 2 a 2 ) RlC, 

(ii) H n h^ (d -2aa 2 )h. 
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Let $ £ OM. be such that Ric($o) > 0- Then the basic relativistic diffusion with diffusivity 
a 2 , started from <3> , explodes with positive probability. 

Note that our choice of constants gives d — 2 a a 2 > and c — 2 a 2 < 0, which gives a 
non-trivial character to conditions (3'). 

4. Conclusion 

We propose in this work a simple probabilistic method to probe certain aspects of the 
singular features of a spacetime, under the form of the existence of incomplete random 
dynamics obtained as random perturbation of the geodesic flow. Theorems [8] and [12] show 
that geodesic completeness can be quite sensitive to random noise under some circumstances. 
This simple approach also highlights some unusual features, compared to the traditional 
studies on geodesic incompleteness. 

- It might be interesting to work on the bundles T : M or OM and not only on the base 
manifold M. 

- The naive arguments given above indicate that causality and initial/boundary con- 
ditions might not be the crux of everything. In that respect, the appearance of a 
dynamic condition on the energy flow happens to be interesting, and needs to be 
compared to the usual pointwise/static energy conditions. 

5. Proofs 

5.1. Proofs of lemma [JJ and proposition [31 

5.1.1. Proof of lemma\J\ • Existence. Lift first arbitrarily the C 2 path (7(5), 7(3)) in T : M 
into a C 1 path T s = (7(5), e(s)) in OM. As -7-7(5) = e (s), there exists some C 1 real-valued 
controls h l ,h 2 ,h 3 and i 1 * 2 , i 1 ' 3 , £ 2 ' 3 defined on [0,T] such that 

3 

r a = H (r s ) + Y l v j (r.)h?(s)+ Mr,)r fc ( s ), 

j=l l^a<6sg3 

where is the canonical vertical vector field on OM generated by the Lie element E a b = 
e« ® e* b ~ e b <g> e* of SO(3) C SO (l,3). Let A be the SO (l, 3)- valued solution of the 
differential equation dA s = — A s E ab £ ab (s). Then the path (^ s )o^ s ^t = (r^sWs^T 

satisfies the conditions of the lemma. 

• Uniqueness. Suppoe (O s )o^s^t is another lift of (7(5), i( s )) 0<S<T to OM satisfying the 
above conditions with some controls g l . As tti(Q s ) = 717 (^s), we need to have O s = ^ S B S , for 
some S0(3) (c S0 (1, 3)) -valued C 2 process (B S ) 0<S<T . Write B s = £ Ka<K3 a ab (s, B s )V ab (B s ), 
identifying here the vector fields V ab to vector fields on SO (3). Then, we have on the one 
hand 

3 

and on the other hand 

3 

e 3 = 4f 3 B a + * 3 B a = H (Q B ) + ^V j (e 3 )ht(8)+ a ab (^B s )v ab (e s ) 

j=l l<a<fes£3 



A PROBABILISTIC VIEW ON SINGULARITIES 



13 



It follows that gi(s) = h^(s), for each j £ {1,2,3}, and all the a ab are identically null, so 
Q s = ^ S B , for some B £ 50(3), and eventually Q s = since O = — To • 

5.1.2. Proof of proposition^ Using the chart F on V provides a trivialization of OV for a 
small enough choice of U. We can thus consider OV as a submanifold of some open set of 
some W. Suppose $ is in OV. Then lemma Q] and Stroock and Varadhan support theorem 
prove that 

P*((* s )(K«sT exits V in F({T} x £/)) > 0. 

To conclude in the general case it suffices to note that one can associate to any pair of points 
(m ,m ), (mi, mi) in T X M a timelike path p : [0, 1] — V M such that (p(0),p(0)) = (m ,m ) 
and (p(l),p(l)) = (mi,mi). 

5.2. Proofs of lemma [9] and theorem [H 

5.2.1. Proof of lemma\B Write (m, (eo,e — 1, 62,63)) for a generic point of OM. As each 
vector field Vj induces no dynamics on M and generates in the SOq(1, 3)-fiber a hyperbolic 
rotation in the 2-plane spanned by e and e^, we have 



Rlc(e tVi $) = -g\R(e , (cosht)e + (sinht)ej) ((cosht)e + (sinht)e^), e 

3 

+ g(^R(e k , (cosht)e + (sinhi)ej) ((cosht)e + (sinht)ej) , e fc 



So we have 



and 



d 2 

V/RlC ($) = — Ric(e^$) = 2 RlC m (e , e ) + 2 RlC m (e i , e j; , 

Ojv t 



(5.1) V/RlC = 4RIC + 4T, 

3=1 

3 _ 
as ^2 RlC m (ej, ej) = 2 T — RlC m (e , e ). The statement of the lemma follows. 

5.2.2. Proof of theorem^ Let $0 be the starting point of the relativistic diffusion in 
we suppose Ric($o) > 0. Under hypothesis (2) we can combine lemmas [7] and [9] and see that 
the diffusion starting from $ explodes with positive probability if there exists a positive 
constant C such that H RlC + 2 a 2 RlC + 2 a 2 T ^ C RlC. As T ^ by hypothesis (1), this 
will be the case if H RlC ^ (C — 2 cr 2 ) RlC; this condition is condition (3). 

5.3. Proof of theorem 1121 The proof consists in checking that we can apply the explosion 
lemma [Til to the non-negative functions / = RlC and h = RlC + U ^ /. As seen in lemma 
[HI the condition Q f ^ cf is equivalent to the inequality 

(5.2) H RlC + 2 a 2 T ^ (c - 2 a 2 ) RlC; 

it follows from condition (3'-i) since T ^ 0. The condition ^ c'h reads 

i? Hy^ 2 Rrc + 2 ff 2 T+^l/ 2 f/^/ l . 

j'=i 3=1 
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By (15. ip and (I3.2p . it is equivalent to 

(5.3) H h + 2 a 2 T < c' h. 

To see that ( 15 .3p follows from condition (3'-ii) notice that the inequality T «C a h is equivalent 
to the inequality RlC+^- ^ U. This condition is implied by condition (2') as we suppose 
R < 0. 
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